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TECHNICAL  NOTE  2834 

FLOW  SURFACES  IN  ROTATING  AXIAL-FLOW  PASSAGES 
By  John  D.  Stanitz  and  Gaylord  0.  Ellis 


SUMMARY 

In  order  to  investigate  the  deviation  of  flow  surfaces  from  their 
assumed  orientation  in  the  usual  type  of  two-dimensional  solution,  three 
dimensional,  incompressible,  nonviscous,  absolute  irrotational  fluid 
motion  is  determined  for  flow  through  rotating  axial-flow  passages 
bounded  by  straight  blades  of  finite  spacing  and  infinite  axial  length 
lying  on  meridional  planes*  Solutions  are  obtained  for  five  passages 
with  varying  blade  spacing  and  hub-tip  ratio-  The  results  are  presented 
in  such  a  manner  as  to  apply  for  all  ratios  of  axial  velocity  to  passage 
tip  speed.  It  is  concluded  that,  for  conditions  in  typical  axial-flow 
blade  rows,  the  deviation  of  flow  surfaces  from  their  assumed  orienta¬ 
tion  in  two-dimensional  solutions  is  small- 


INTRODUCTION 

A  flow  surface  in  the  passage  between  two  blades  of  a  compressor 
or  turbine  is  generated  by  the  motion  through  the  passage  of  any  fluid 
line  consisting  of  the  same  fluid  particles  and  extending  from  one 
boundary  to  another  in  a  plane  normal  to  the  axis  of  rotation-  In  two- 
dimensional  analyses  of  flow  in  compressors  and  turbines,  the  fluid 
motion  is  usually  assumed  to  occur  on  flow  surfaces  that  are:  (l)  sur¬ 
faces  of  revolution  about  the  axis  of  the  turbomachine  (blade -to -blade 
solutions,  references  1  and  2,  for  example)  or  (2)  mean  passage  surfaces 
that  are  congruent  with  the  mean  blade  surfaces  (hub-to-shroud  solutions 
references  3  and  4,  for  example) -  Actually,  the  flow  surfaces  deviate 
from  the  orientation  assumed  for  the  two-dimensional  solutions  and,  in 
the  direction  of  flow,  become  progressively  more  tilted  and  distorted. 
This  deviation  of  the  flow  surfaces  from  their  assumed  Orientation  is 
caused  by  spanwise  variations  of  blade  loading  and,  in  rotating  blade 
rows,  by  rotation  of  the  fluid  particles  relative  to  the  passage  in  a 
plane  normal  to  the  axis  of  the  blade  row-  This  rotation  is  required 
to  maintain  the  rotational  or  irrotational  character  of  the  absolute 
fluid  motion. 
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The  deviation  of  flow  surfaces  is  considered  in  reference  5,  hut  no 
attempt  is  made  to  estimate  the  magnitude  of  this  phenomenon. .  An  ana¬ 
lytical  investigation  has  therefore  been  made  at  the  NACA  Lewis  lab ora- 
tory  in  order  to  determine  the  magnitude  of  this  deviation  in  rotating 
axial-flow  passages.  The  axial-flow  passages  in  this  investigation  are 
hounded  hy  straight  blades  of  finite  spacing  and  infinite  axial  leng 
lying  on  meridional  (axial-radial)  planes.  The  solutions  have  been  made 
for  three-dimensional,  incompressible,  nonviscous,  absolute  irrotational 
fluid  motion  over  a  range  of  blade  spacings  and  hub -tip  ratios.  These 
solutions  do  not  investigate  the  effect  of  spanwise  distribution  of  blade 
loading,  which  was  considered  of  secondary  importance.  (Note  that,  as  m 
rectangular  elbows  with  potential  flow,  uniform  spanwise  loading  has  no 
effect  on  the  deviation  of  flow  surfaces.)  Likewise,  the  effects  of  - 
compressibility  have  not  been  investigated  because,  as  clearly  indicate 
by  the  correlation  equations  in  reference  1,  the  eddy  flow,  which  causes 
the  flow  surfaces  to  deviate,  is  little  affected  by  compressibility. 

The  results  are  presented  in  such  a  manner  as  to  apply  to  any  ratio  of 
blade-tip  speed  to  axial  velocity  of  the  fluid. 


METHOD  OF  SOLUTION 

The  method  of  solution,  including  the  relaxation  solution  of  the 
differential  equation  of  flow  and  the  superposition  of  solutions,  is 
developed  in  this  section. 


Preliminary  Considerations 


Assumptions.  -  The  absolute  flow  is  assumed  to  be  irrotational- 
The  fluid  is  assumed  to  be  nonviscous  and  incompressible.  The  fluid 
motion  is  three  dimensional  and  is  steady  relative  to  the  rotating 
passage. 


Coordinate  system  and  velocity  components .  -  The  cylindrical 
coordinate  system  r,  0,  z  is  shown  in  figure  1.  (All  symbols  are 
defined  in  appendix  A) .  The  linear  coordinates  r  and  z  are  expressed 
as  ratios  of  the  blade -tip  radius.  Thus,  for  example,  the  radius  r 
at  the  blade  tip  is  unity.  The  coordinate  system  is  fixed  relative  o 
the  passage  which  rotates  about  the  z-axis  in  the  positive  direction 
according  to  the  right-hand  rule. 


The  velocity  components  u,  v,  w  relative  to  the  coordinate  system 
in  the  r,  0,  z  directions,  respectively,  are  also  shown  in  figure  1 ■ 
The  velocity  components  and  the  blade  speed  are  expressed  as  ratios  of 
the  blade -tip  speed.  Thus,  for  example,  the  blade  speed  at  any  radius 
is  equal  to  r  and  the  absolute  tangential  velocity  component  becomes 

(v  +  r)  . 
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Type  of  passage  geometry.  -  The  rotating  axial-flow  passages  in 
this  investigation  are  infinitely  long.  Each  passage  is  hounded  hy  a 
huh  and  casing  of  constant  radius,  respectively,  and  hy  straight  blades 
of  finite  spacing  and  infinite  length  lying  on  meridional  (axial-radial) 
planes.  The  hlade  inlet  is  considered  to  he  at  minus  infinity  in  the 
z-direction  and  the  hlade  exit  at  plus  infinity.  Under  these  circum¬ 
stances  the  flow  is  uniform  in  the  z-direction  at  the  region  investi¬ 
gated  (near  the  origin,  z  =  0)  and  the  hlade  loading  is  zero.  Thus, 
effects  of  hlade  loading  on  deviation  of  the  flow  surfaces  are  not  inves¬ 
tigated  in  this  report.  These  effects  are  considered  of  secondary 
importance . 

Superposition  of  solutions.  -  For  the  passage  geometry  just 
described,  the  incompressible  flow  solution  can  he  separated  into  two 
parts:  (l)  the  rotating  or  eddy-flow  solution  in  the  rotating  passage 
with  no  through  flow  and  (2)  the  through-flow  solution  in  the  station¬ 
ary  passage  with  no  eddy  flow.  The  eddy -flow  solution  does  not  change 
in  the  z-direction  and  is  therefore  two  dimensional.  The  through-flow 
solution  is  a  uniform  axial  velocity  w.  Various  percentages  of  the 
two  solutions  can  he  combined  hy  linear  superposition  to  obtain  new 
solutions  for  different  ratios  of  axial  velocity  to  blade-tip  speed, 
that  is,  for  different  values  of  w. 


Eddy -Flow  Solution 

The  eddy -flow  solution  is  two  dimensional  and  lies  in  the  r 6 -plane. 

Continuity.  -  A  fluid  particle  on  the  r6-plane  is  shown  in  figure  2 
From  continuity  considerations 

'  (1) 

A  stream  function  if  satisfies  equation  (l)  if  defined  as 

ru  (2a) 

g  =  -v  (2b) 

Irrotational  absolute  motion.  -  For  irrotational  absolute  motion, 
the  circulation  of  the  absolute  velocity  around  the  fluid  particle  in 
figure  2  is  zero,  and  therefore 


Cr  +  v)r]  -  0 
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or 


dv  v  1  du  , 

Br  +  r  "  r  3©  - 


(3a) 


■which..?  after  substitution  of  equation  (2),  becomes 

d2^  1  cty  ,  _1_  d2ty  _  2 

Sr2  r  ^  r2  be2 


(3b) 


Equation  (3b)  is  the  differential  equation  of  flow  that  determines  the 
distribution  of  t  for  the  eddy-flow  solution  in  the  r 6 -plane. 


Transformation  of  coordinates.  -  In  order  to  solve  equation  (3)  by 
relaxation  methods,  it  is  convenient  to  transform  the  r0-plane  onto  the 
£ 0-plane  by  means  of 

|  =  In  r  (4=) 


from  which  equation  (3b)  becomes 


d2i|r 

bl2  be2 


(5) 


Relaxation  solution.  -  Equation  (5)  is  solved  by  relaxation  methods 
(references  6  and  7,  for  example)  to  satisfy  the  specified  boundary 
conditions.  For  the  eddy-flow  solutions,  there  is  no  flow  through  the 
passage  so  that  i  is  zero  along  the  hub,  shroud,  and  blade  surfaces. 

In  the  -plane  these  boundaries  form  a  rectangle  within  which  is  placed 
a  grid  of  equally  spaced  points .  At  each  of  these  grid  points  the  value 
of  ilr  required  to  satisfy  equation  (5)  in  finite  difference  form  is 
determined  by  relaxation  methods.  The  size  of  the  grid  spacing  varies 
among  examples  and  will  be  indicated  later.  The  values  of  *  at  the 
grid  points  were  relaxed  to  a  unit  change  in  the  fifth  decimal.  The. 
velocity  components  are  obtained  from  the  distribution  of  ijr  according 
to  equation  (2).  The  streamlines  of  the  eddy  flow  in  the  r@-plane  are 
lines  of  constant  \|r. 


Combined  Solutions 

For  the  eddy-flow  solutions  on  the  r0-plane,  the  fluid  rotates 
relative  to  the  passage  walls  in  a  direction  opposite  to  that  of  the 
blade  rotation.  This  fluid  motion  is  the  same  for  all  planes  normal  to 
the  z-axis.  For  the  through-flow  solution  the  axial  velocity  w  is 
everywhere  constant.  These  two  linear  solutions  can  be  superposed  to 
obtain  solutions  for  three-dimensional  flow  through  rotating  axial-flow 
passages. 
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It  is  desired  to  determine  the  flow  surface  generated  by  the 
motion  of  any  fluid  line  that  extends  between  boundaries  in  the  r 6 -plane 
and  always  consists  of  the  same  fluid  particles .  This  fluid  line 
rotates  with  the  fluid  in  the  r0 -plane  and  the  surface  that  it  generates 
depends  on  the  velocity  w  with  which  it  moves  in  the  axial  direction 
through  the  passage •  Examples  of  such  flow  surfaces  are  shown  in  fig¬ 
ure  3.  The  shape  of  these  surfaces  can  be  indicated  on  the  r9-plane 
alone  by  plots  of  the  intersections  of  these  surfaces  with  the  r 9 -plane 
at  equal  increments  of  z.  These  intersections  are  the  positions  of  the 
fluid  lines  on  the  r0-planes  at  these  values  of  z.  If,  instead  of 
increments  of  z,  however,  fluid  lines  are  plotted  on  the  r 9 -plane  for 
increments  of  the  absolute  angle  a  that  the  passage  has  rotated  about 
the  z-axis,  these  fluid- line  positions  apply  for  all  values  of  w.  For 
a  given  value  of  a,  the  value  of  z  then  depends  on  w  and  this  rela¬ 
tion  is  given  by 

z  =  aw  (6) 

Thus  the  results  of  the  combined  solutions  are  plotted  as  fluid-line 
positions  in  the  r@ -plane  for  equal  increments  of  a  and  these  results 
(fig.  4(a),  for  example)  apply  for  all  ratios  w  of  through-flow  velocity 
to  blade-tip  speed.  The  three-dimensional  flow  surfaces  in  figure  3 
correspond  to  the  fluid-line  positions  shown  in  figure  4(a)  for  w 
equal  to  0.6. 


NUMERICAL  EXAMPLES 


The  results  for  three-dimensional  flow  through  five  rotating  axial- 
flow  passages  are  presented  in  figures'  4  to  8.  The  results  are  pre¬ 
sented  in  the  r0 -plane  by  streamlines  \|r*  of  the  eddy-flow  solution  and 
by  fluid-line  positions  a  of  the  three-dimensional  flow  surfaces.  The 
streamlines  are  designated  by  T|r*,  which  is  defined  as 


**  =  — i- 

^min 


(7) 


where  the  subscript  min  refers  to  the  algebraic  minimum  value  of  f 
so  that  V  varies  from  zero  along  the  boundaries  to  1*0  at  the  point  of 
minimum  i|r.  The  fluid-line  positions  in  the  r9-plane  are  indicated  for 
various  values  of  the'  absolute  angle  a  that  the  passage  has  rotated 
about  the  z-axis  from  its  initial  position  (a  =  0)  at  which  the  fluid¬ 
line  positions  are  radial  or  circumferential  lines .  In  conformity  with 
reference  5,  flow  surfaces  with  initial  fluid- line  positions  that  are 
circumferential  or  radial  lines  are  designated  S-|_-  or  S2-surfaces, 
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Passage  configurations*  -  The  geometry  of  the  five  axial-flow 
passages  investigated  is  described  in  table  I.  • 


TABLE  I  -  GEOMETRY  OF  AXIAL-FLOW  PASSAGES 


Example 

Hub -tip 
ratio 

Blade  spacing 
A0 

Grid 

spacing 

rh 

radians 

deg 

l( standard) 

0.70000 

0.17834 

a 

0 

s 

(A0/8)  =  0.02229 

II 

.70000 

.08917 

5°  7' 

(A0/8)  =  .01115 

hi 

. 70000 

.35667 

20°26' 

(A0/8)  =  .04458 

IV 

•50105 

. 17834 

10° 13 ' 

(A0/8)  =  .02229 

V 

. 89453 

. 17834 

10° 13' 

(A0/l6)=  .01115 

The  results  of  the  standard  solution,  example  I,  are  compared  with  the 
results  of  examples  II  and  III  to  determine  the  effect  of  varying  the 
blade  spacing  A 0  with  constant  hub -tip  ratio  r^.  The  results  of 
example  I  are  also  compared  with  examples  IV  and  V  to  determine  the. 
effect  of  varying  r^  with  A0  constant.  The  grid  spacings  used  in 
the  relaxation  solutions  are  given  in  the  last  column  of  the  table- 

Standard  solution.  -  Results  for  the  standard  solution  (example  i) 
are  presented  in  figure  4.  In  figure  4(a)  are  shown  fluid-line  posi¬ 
tions  of  the  central  flow  surfaces  for  various  values  of  the  angle  a. 

The  central  flow  surfaces  are  defined  as  those  surfaces  for  which  the 
fluid  lines  pass  through  the  point  of  minimum  i|r,  that  is,  f  =  1.0. 

At  this  point,  values  of  u  and  v  are  both  zero  so  that  the  central 
flow  surfaces  pivot  about  a  straight  line  in  the  z-direction  through 
this  point. 

In  figure  4(b)  are  shown  the  fluid-line  positions  of  off-center 
Sp-surfaces  for  various  values  of  the  angle  a.  For  any  off-center  flow 
surface,  the  envelope  of  the  fluid-line  positions  for  various  values  of 
a  is  a  streamline.  This  fact  is  clearly  shown  by  the  upper  Sp-surface 
in  figure  4(b)  which  is  tangent  to  the  streamline  0.8. 

Fluid-line  positions  of  off-center  S2-surfaces  for  various  values 
of  a  are  shown  in  figure  4(c).  Finally,  in  figure  4(d),  are  shown 
fluid-line  positions  of  the  central  Sp-surface  for  a  wide  range  of  a. 

As  a  increases,  the  surface  becomes  progressively  more  distorted  because 
its  velocity  in  the  r 0 -p lane  along  the  boundaries  near  the  corners  is 
low,  becoming  equal  to  zero  at  the  corners,  whereas  the  velocities  along 
most  of  the  other  eddy-flow  streamlines  approach  a  wheel-type  distribu¬ 
tion  with  zero  velocity  at  —  1.0.  It  is  concluded  that,  for  large 
values  of  the  absolute  angle  a,  the  flow  surfaces  become  greatly 
distorted. 
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Solutions  for  effect  of  Id  lade  spacing.  -  Examples  II  and  III  are 
presented  in  figures  5  and  6.  These  figures,  together  with  figure  4(a), 
indicate  the  shapes  of  the  central  flow  surfaces  for  three  "blade  spacings 
A0  with  the  same  hub -tip  ratio  r^.  The  general  appearance  of  the 
central  S^-surfaces  is  similar  for  examples  I  and  II,  and  in  example  III 
the  Sq-  and  S2-surfaces  are  similar.  Reasons  for  these  similarities  are 
given  in  DISCUSSION  OF  RESULTS. 

Solutions  for  effect  of  hub -tip  ratio.  -  Examples  IV'  and  V  are 
presented  in  figures  7  and  8.  These  figures,  together  with  figure  4(a), 
indicate  the  shapes  of  the  central  flow  surfaces  for  three  hub -tip  ratios 
rh  with  the  same  blade  spacing  AS.  The  general  appearance  of  the 
central  Sq-surfaces  of  examples  I  and  IV  is  similar,  and  the  Sq-  and 
Sg-surfaces  of  example  V  are  similar  to  the  Sg-  and  Sq-surfaces,  respec¬ 
tively,  of  example  I.  Also,  it  is  noted  that  the  central  Sq-  and  S2- 
surfaces  of  examples  II  and  IV  are  similar  in  general  appearance.  Rea¬ 
sons  for  these  similarities  are  given  in  DISCUSSION  OF  RESULTS- 


DISCUSSION  OF  RESULTS 

Some  of  the  results  presented  in  figures  4  to  8  are  discussed,  and 
the  deviations  of  the  flow  surfaces  from  their  initial  positions  for  a 
'equal  to  zero  in  the  rS-plane  are  investigated. 

,  Typical  value  for  a.  -  The  results  in  figures  4  to  8  are  presented 
as  fluid-line  positions  in  the  rS-plane  for  even  increments  of  a..  As 
already  defined,  a  is  the  absolute  angle  that  the  axial-flow  passage  , 
has  rotated  about  the  z-axis,  with  a  equal  to  zero  when  the  initial 
position  of  the  fluid  line  is  a  circumferential  line  (Sq-surface)  or 
radial  line  (S2 -surface)  in  the  rS-plane.  This  angle  a  is  related  to 
the  geometry  and  operating  conditions  of  the  axial-flow  passage  by 
equation  (6).  In  order  to  determine  a  typical  value  for  a,  an  axial- 
flow  stage  is  considered  with 

w’  =  550  ft/sec 

z’  =  0.12  ft  (1-44  in.) 

0)' '  =  838  radians/sec  (8000  rpm) 

where  the  prime  superscript  indicates  dimensional  quantities.  Equa¬ 
tion  (6)  becomes 
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so  that  a  typical  value  for  a  is  approximately  0.2.  From  figures  4  to 
8  it  is  therefore  concluded  that  the  deviation  of  flow  surfaces  in 
typical  axial-flow  blade  rows  is  not  large.  This  conclusion  is  further 
strengthened  if  the  fluid-line  position  for  <x  equal  to  zero  is  con¬ 
sidered  to  occur  halfway  through  the  blade  row.  Then  the  maximum 
deviation  of  the  surface  from  its  position  at  a  equal  to  zero  is 
reduced  by  approximately  one  half. 

For  blade  rows  (not  necessarily  axial  flow)  with  relatively  large 
dimensions  in  the  direction  of  flow,  such  as  radial-  and  mixed-flow 
impellers,  the  deviations  of  the  flow  surfaces  must  be  large.  However, 
even  these  large  deviations  do  not  invalidate  the  two-dimensional  solu¬ 
tions  completely,  because,  as  shown  in  reference  8,  at  many  positions 
in  the  passage  the  velocity  components  of  major  importance  are  much  the 
same  for  two-  and  three-dimensional  solutions. 

Deviation  of  flow  surfaces .  -  The  deviation  of  flow  surfaces  from 
their  initial  orientation,  given  by  fluid-line  positions  in  the  rS-plane 
at  a  equal  zero,  can  be  described  by  three  factors  (fig*  9) :  (!) 

displacement,  in  the  r@ -plane ,  of  the  tangent  point  between  the  fluid 
line  and  the  tangent  streamline;  (2)  rotation,  in  the  r0-plane,  of  the 
fluid  line  about  this  tangent  point;  and  (3)  distortion  or  bending  of 
the  fluid  line  in  the  r0 -plane .  The  displacement  of  the  tangent  point 
is  determined  by  its  motion  in  the  r@-plane  along  the  streamline  with 
which  the  fluid  line  is  tangent.  This  displacement  for  off-center  Sp- 
and  S2 -surfaces  is  indicated  in  figures  4(b),  4(c),  and  9  and  will  not 
be  discussed  further.  For  central  flow  surfaces,  which  will  be  con¬ 
sidered  exclusively  hereinafter,  the  tangent  point  (center  point)  does 
not  move  and  the  displacement  is  zero. 

The  rotation  of  central  flow  surfaces  will  be  measured  by  the  angle 
P  -  P0  which  the  tangent  to  the  fluid  line  at  its  center  point  rotates 
in  the  r0 -plane  from  its  initial  position  Pq  at  a  equals  zero 
(fig.  9) .  The  angle  p  is  measured  clockwise  from  the  radial  direction 
so  that  Po  is  90°  for  Sp-surf aces  and  0°  for  Sg-surfaces. 


The  distortion  of  the  flow  surfaces  will  be  discussed  qualitatively. 


Rotation  of  flow  surfaces.  -  The  rotation  of  central  flow  surfaces 
is  measured  by  the  angle  p  -  Pq  introduced  in  the  preceding  section. 
This  angle  can  be  measured  in  figures  4(a)  and  5  to  8.  However,  an 
equation  has  been  developed  (appendix  B)  by  which  the  rotation  P  -  Po 
can  be  determined  directly  from  a  and  a  parameter  A,  which  is '  the 

value  of  — 5-  at  the  center  point. 

dr2 
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For  central  Sp-surfaces 

tan(p  -  PQ)1  =  -  tan  [a  a/A(2  -  A)]  (7a) 

and  for  central  S2 -surfaces 

tan(p  -  P0)2  -  VS  —  tan  j^a  /\/a(2  -  A)J  (7b) 

so  that  the  rotations  of  the  two  types  of  central  flow  surface  are 
related  by 

tan(p  -  P0) i  =  2'”)  tan^P  "  P0^2  (7c) 

In  particular,  for  A  equal  to  1.0, 

(P  -  P0)l  "  (P  “  Po)2  “  a  (7d) 


so  that  the  rotation  of  both  flow  surfaces  are  equal  to  the  rotation  a 
of  the  passage  about  the  z-axis. 


As  will  be  discussed  later  in  this  section,  the  parameter  A  is 


In  r 


primarily  a  function  of  — 259”*  ^ 


In  rh 


1.0  and  A8  is  finite,  v  and  are  zero  at  the  center  point  and 

equation  (3a)  gives 

dv  _  2 


Ae 


is  zero,  that  is,  if  ri 


is 


from  which 


Also,  if 


In  r. 

~ A 6 

dv 


is  infinite,  that  is,  if  A0  is  zero  and  r^  is  less 


than  1.0,  is  zero  at  the  center  point  so  that 


A  =  11=  -  |l 

Sr2  Sr 


0 


in  rh 

Thus,  the  parameter  A  varies  between  0  and  2.0  as  — varies 
between  -oo  and  0.  For  A  equal  to  0,  equations  (7a)  and  (7b)  become 
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tan(p  -  Pq^I  =  “2<x  ^ & ^ 

and  v 

tan(p  -  Po^2  “  0 

Likewise,  for  A  equal  to  2.0,  equations  (7a)  and  (7b)  become 

tan(p  -  Pq) i  =0  (?g) 

and 

tan(P  -  Pq)2  =  -2<x 


The  rotation  (p  -  p0)x  of  central  S^-surfaces  has  been  computed 

by  equations  (7a),  (7e),.and  (7g)  and  is  plotted. in  figure  10  as  a  func¬ 
tion  of  a  for  various  values  of  A.  The  rotation  (p  -  Po'2  °* 
central  S2-surfaces  is  also  given  by  figure  10  if  the  curves  of  constant 
A  are  numbered  in  reverse  order-  Thus,  discussions  relating  to  the 
rotation  of  central  Sq- surfaces  with  parameter  A  equal  to  x  also 
apply  to  the  rotation  of  central  Sg-surfaces  with  A  equal  to  (2  -  x) . 

In  figure  10  the  curve  for  A  equal  to  zero  is  asymptotic  to  rt/2, 

or  1.5708.  For  this  value  of  A,  the  passage  width  is  zero  {A6  =  0) , 
and  the  central  S]_-surface  cannot  rotate  more  than  rt/2  radians.  For 
A  equal  to  2.0,  the  rotation  (p  -  Po)i  is  zero  at  a11  values  of  a’ 

For  this  value  of  A,  the  passage  height  is  zero  (r^  =  1*0),  and  the 

central  S]_-surface  cannot  rotate.  As  indicated  by  equation  (7d),  a 
linear  relation  exists  between  (P  -  Pq)i  an<^  a  ^ or  ^  e<lual  "^°  1*0* 
As  will  be  shown  later  in  this  section,  for  this  value  of  A  the  average 
passage  width  is  approximately  equal  to  the  passage  height  (example  III, 
fig.  6),  and  both  the  central  Sx-  and  S2-surfaces  rotate  at  the  same 
rate  as  the  passage  itself,  but  in  the  opposite  direction.  For  the 
remaining  values  of  A,  the  curves  in  figure  10  have  inflection  points 
at  (P  -  Pq)!  equal  to  it/2,  rt,  and  so  forth-  For  values  of  A  less 

than  1.0,  the  rate  of  change  of  (P  -  Pq)i  with  a  is  minimum  at 
(p  -  P0)-l  equal  to  |,  3|,  and  so  forth,  and  is  maximum' at  (p  -  P0)i 
equal  to  *,  2it,  and  so  forth.  For  values  of  A  greater  than  1.0,  the 
reverse  is  true.  In  all  cases,  the  rate  of  change  of  (P  -  P0)i  Wlth 
a  is  greatest  when  the  tangent  to  the  fluid  line  at  its  center  point  is 
oriented  in  the  direction  of  minimum  distance  between  passage  walls  and 
is  least  when  the  tangent  is  oriented  normal  to  the  direction  of  minimum 
distance.  This  observation  is  reasonable  because,  as  indicated  by  the 
streamline  spacing  for  examples  I  to  V,  the  gradient  of  the  velocity 
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component  normal  to  the  tangent  of  the  fluid  line  at  its  center  point, 
which  velocity  gradient  causes  the  fluid  line  to  rotate  about  its  center 
point,  is  maximum  when  the  tangent  is  oriented  in  the  direction  of  min¬ 
imum  distance  between  passage  walls  and  is  minimum  when  the  tangent  is 
oriented  normal  to  the  direction  of  minimum  distance. 

The  parameter  A,  which  determines  the  rotation  of  the  central  flow 

In  rh 

surfaces,  is  primarily  a  function  of  the  ratio  — values  A 
for  examples  I  to  V  have  been  obtained  from  the  relaxation  solutions  and 
are  given  in  table  XX  together  with  the  values  of 


TABLE  II  -  VALUES  OF  PARAMETER  A 


Example 

A 

In  r^ 

H. 

ag 

I 

0.314 

-2 

ii 

.134 

-4 

hi 

1.020 

-i 

IV 

.240 

.5! 

o 

V 

1.535 

5 

“8 

—  rk.  These  values  of  A  and  are  plotted  in  figure  11.  As 

a  e  ag  t 

previously  discussed,  the  parameter  A  is  equal  to  2  and  zero  for 

^J1  rk  equal  to  zero  and  -<»,  respectively.  It  can  be  shown  analytically 
A G  In  rh 

that  the  curve  in  figure  11  has  zero  slope  for  —?~z —  equal  to  zero. 

In  rh 

As  — — =  varies  from  zero  to  -o »,  the  passage  geometry  in  the  r 6 -plane 

varies  from  a  wide  shape  with  zero  height  in  "the  r-direction  to  a  tall 

In  r^ 

shape  with  zero  width  in  the  ©-direction.  For  — ~/\§~  equal  to  -1.0, 

the  passage  geometry  is  square  in  the  | ©-plane  and  the  average  passage 
width  in  the  r ©-plane  is  approximately  equal  to  the  passage  height 
(example  III,  fig.  6). 

In  rh 

The  parameter  A  is  a  function  of  — ~^q~  because  passages  with 
the  same  value  of  — have  geometrically  similar  boundaries  in  the 

/  \  m  I  m 

£0-plane,  where  the  solution  of  equation  (5) ,  which  solution  determines 

A,  is  obtained.  The  right  side  of  equation  (5)  indicates,  however,  that 

the  solution  of  equation  (5) ,  and  therefore  the  value  of  A,  depends  not 

.  in  r^ 

only  on  the  passage  shape  in  the  £© -plane,  that  is  on  — &q~>  also 

on  the  corresponding  values  of  r  at  each  value  of  Thus,  the 
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value  of  A  must  also  depend  on  the  hub-tip  radius  ratio  rh-  However, 
figure  11  shows  that  rh  has  only  a  small  effect  on  A  for  the  range 
of  rh  investigated. 

For  the  range  of  a  investigated  by  the  numerical  examples 
(0  <  a  <  1.0),  figure  10  shows  that  the  variation  in  (p  -  P0)i  with  a 
is  similar  for  0  <  A<  0.3  (also  compare  examples  I,  II,  and  IV  in 
figs.  4(a),  5,  and  7,  respectively);  and,  if  the  curves  of  constant  A 
are  numbered  in  reverse  order,  figure  10  indicates  the  variation  in 
p0)g  with  a  is  similar  for  1.7  <  A  <  2.0.  In  both  cases  the 

rotation  of  the  central  flow  surfaces  is  similar  for  the  specified  range 

.  '  -  ■  ,  ln  rh 

of  A  because  for  this  range  the  corresponding  values  of 

(fig.  ll)  are  such  that  the  passage  walls  parallel  to  the  initial  posi¬ 
tions  (a  =  0)  of  the  central  flow  surfaces  are  too  far  removed  to 
exert  an  important  influence  on  the  rotation  (P  -  P0) ,  which  is  therefore 
affected  primarily  by  the  angle  a.  It  is  therefore  concluded  that: 

(1)  For  values  of  —ggT  algebraically  less  than  -2,  the  rotation  of 
central  S^-surfaces  is  about  the  same  for  a  less  than  1.0;  and  (2)  for 

values  of  - algebraically  greater  than  -0.5,  the  rotation  of  cen- 

A  6 

tral  S2-surfaces  is  about  the  same  for  a  less  than  1.0. 

Distortion  of ! flow  surfaces .  -  Factors  affecting  the  distortion  of 
the  flow  surfaces  are  evident  from  figures  4  to  8.  In  general,  a  sur¬ 
face  becomes  distorted  if  (1)'  the  fluid  line  that  generates  the  surface 
approaches  the  vicinity  of  a  corner  in  the  r 6-plane  and  (2)  the  center, 
or  tangency  point  of  the  fluid  line,  moves  closer  to  one  of  the  passage 
boundaries.  The  relative  importance  of  these  factors  depends  on  the 
particular  passage  geometry  and  the  orientation  of  the  flow  surface* 

From  figure  4(b) ,  if  the  fluid  line  of  the  off-center  flow  surfaces  is 
initially  oriented  (a  =  0)  normal  to  the  longer  side  of  the  passage 
boundary,  the  first  factor  is  of  major  importance.  From  figure  4(c),  if 
the  fluid  line  of  the  off-center  flow  surfaces  is  initially  oriented 
parallel  to  the  longer  side,  the  second  factor  is  most  important.  For 
central  flow  surfaces,  only  the  first  factor  exists. 

It  is  clearly  evident  from  figure  4(d)  that  for  values  of  a  con¬ 
siderably  larger  than  1-0  the  flow  surfaces  become  greatly  distorted. 

Effect  of  —At.  _  Table  II  and  figure  11  indicate  an  approximate 

_  Ae  In  rh  • 

correlation  of  the  parameter  A  with  the  ratio  2\6~‘  Because  A  is 

an  important  parameter  in  the  calculation  of  the  rotation  (P  -  Pq)  by 
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In  rh 

equation  ( 7 ),  this  correlation  suggests  that  — is  an  important 
parameter  affecting  the  shape  of  the  flow  surfaces. 


In  table  II  the  value  of 


In  r^ 


is  nearly  the  same  for  examples  II 


and  IV,  and  a  comparison  of  these  examples  in  figures  5  and  7  indicates 
great  similarity  in  the  shape  of  the  flow  surfaces.  Also,  in  table  II, 

lll-r—  for  example  I  and  ,  —  ■  for  example  V  are  of  the  same  general 
A9  In  r^ 

order  of  magnitude.  A  comparison  of  the  Sq-  and  S2 -surfaces  of  example  I 
in  figure  4(a)  with  the  S2-  and  S^-surfaces ,  respectively,  of  example  V 

in  figure  8  indicates  considerable  similarity.  It  is  therefore  concluded 

In  rh 

that,  for  the  same  value  of  — &q~>  or  inverse>  axial -flow  passages 
of  the  type  investigated  have  similar  shapes  of  flow  surfaces. 


SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

Three-dimensional,  incompressible,  noiiviscous,  absolute  irrotational 
fluid  motion  is  investigated  for  flow  through  rotating  axial-flow  pas¬ 
sages  bounded  by  straight  blades  of  finite  spacing  and  infinite  axial 
length  lying  on  meridional  planes.  Solutions  are  obtained  for  five  pas¬ 
sage  geometries  described  by  various  ratios  of  the  logarithm  of  the  hub- 

ln  rh 

tip  ratio  divided  by  the  blade  spacing  — ,  and  the  results  are 
presented  in»  such  a  manner  as  to  apply  for  all  ratios  of  axial  velocity 
to  passage  tip  speed. 

The  solutions  are  used  to  determine  the  deviation  of  flow  surfaces 
from  their  assumed  orientation  in  the  usual  type  of  two-dimensional 
solution.  This  deviation  is  shown  by  the  fluid-line  positions  (inter¬ 
sections  of  the  flow  surfaces  with  the  r 6 -plane)  for  equal  increments  of 
the  angle  a  that  the  passage  rotates  about  the  z-axis  as  the  flow 
surface  deviates  from  its  initial  orientation.  The  deviation  is  con¬ 
sidered  to  consist  in  (l)  displacement  in  the  r0-plane  of  the  center 
point  of  the  fluid  line,  (2)  rotation  in  the  r 6-plane  of  the  fluid  line 
about  its  center  point,  and  (3)  distortion  of  the  fluid  line  in  the 
r0 -plane. 

Two  types  of  flow  surface  are  considered:  Sq-  and  Sg -surfaces 
initially  oriented  along  circumferential  and  radial  lines,  respectively, 
in  the  r6-plane.  The  surfaces  are  central  flow  surfaces  if  they  pass 
through  the  point  of  zero  relative  velocity  at  the  passage  center  in  the 
re-plane;  otherwise,  the  surfaces  are  off-center  flow  surfaces. 


14 


NACA  TN  2834 


Some  results  of  the  numerical  examples  are: 

]_.  The  central  flow  surfaces  rotate  relative  to  the  passage  about 
a  straight  axial  line  through  the  point  of  minimum  stream  function  near 
the  center  of  the  passage  in  the  r 6 -plane. 

2.  For  any  off-center  flow  surface  the  envelope  of  the  fluid-line 
positions  in  the  r0-plane  for  various  values  of  a  is  a  streamline. 

Some  conclusions  resulting  from  the  numerical  examples  are: 

1.  For  values  of  a.  corresponding  to  conditions  in  typical  axial- 
flow  blade  rows,  the  deviation  of  flow  surfaces  is  not  large. 

2.  For  values  of  a  less  than  1.0  radian,  the  rotation  of  the 

.•I  ^h 

central  S^- surface  is  about  the  same  in  all  passages  for  which 
is  algebraically  less  than  -2.0. 

3.  For  values  of  a.  less  than  1.0  radian,  the  rotation  of  the 

central  Sg -surface  is  about  the  same  in  all  passages  for  which 
is  algebraically  greater  than  -0.5. 

4.  In  general,  a  flow  surface  becomes  distorted  if  (a)  the  fluid 
line  that  generates  the  surface  approaches  the  vicinity  of  a  corner  in 
the  r@-plB.ne  and  (b)  the  center  point  of  the  fluid  line  moves  closer  to 
one  of  the  passage  boundaries. 

5.  For  values  of  a  considerably  greater  than  1.0  radian,  the  flow 
surfaces  become  greatly  distorted. 

6.  For  the  same  value  of  — axial-flow  passages  of  the  type 
investigated  have  similar  shapes  of  flow  surfaces. 

Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  August  1,  1952 


2617 


2617 


MCA  TN  2834 


15 


APPENDIX  A 


SYMBOLS 

The  following  symbols  are  used  in  this  report.  All  symbols  are 
dimensionless,  unless  otherwise  specified.  Velocities  are  expressed  as 
ratios  of  the  passage-tip  speed;  distances  are  expressed  as  ratios  of 
the  passage-tip  radius. 

A  parameter,  which  is  the  value  of 

(\jr*  s  l.o)  about  which  the  central  flow  surfaces  rotate, 
equation  (B3a)  of  appendix  B 

r,0,z  cylindrical  coordinates  relative  to  rotating  passage  (fig.  l) 

S  flow  surface  generated  by  motion  through  passage  of  any  fluid 

line  consisting  of  the  same  fluid  particles  and  extending 
from  one  boundary  to  another  in  r 6 -plane 

s  arc  length  along  flow  surface  in  r 6 -plane 

u, v,w  relative  velocity  components  in  r,  6,  z  directions,  respectively, 
(fig-  1) 

a  absolute  angle  that  passage  has  rotated  about  z-axis  from 

initial  position  at  -which  fluid  lines  for  and  Sg-surfaces 
are  circumferential  and  radial  lines,  respectively,  in 
r Q -plane 

p  angle  of  tangent  to  fluid  line  at  its  center  point  in  r 0-plane, 

measured  clockwise  from  radial  direction 

A 6  blade  spacing  in  r0-plane 

|  transformed  coordinate,  equation  (4) 

stream  function  in  r0-plane,  equation  (2) 

r*  stream  function  \|r  divided  by  equation  (7) 

Q  relative  angular  velocity  of  elemental  arc  ds,  of  central  flow 

surface,  rotating  about  point  =  1*0  in  r0-plane,  expressed 
as  ratio  of  cd1 

absolute  angular  velocity  of  passage  about  z-axis,  dimensional 
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Subscripts : 

h  hub  (so  that  rh  is  hub -tip  ratio) 
min  minimum 

0  initial  position,  when  a  equals  zero 

1  flow  surface  with  circumferential  line  for  initial  position  of 

fluid  line  in  r0 -plane 

2  flow  surface  with  radial  line  for  initial  position  of  fluid  line 

in  r0-plane 

Superscript: 


f 


dimensional  quantities 
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APPENDIX  B 


1* 


ROTATION  OF  CENTRAL  Sx-  AND  S2 -SURFACES  ABOUT  THEIR  CENTER  POINT 

IN  re-PLANE 


If  Q  is  the  relative  angular  velocity,  expressed  as  a  ratio  of 
the  absolute  angular  velocity  of  the  passage  about  the  z-axis,  of  an 
elemental  arc  length  ds  rotating  about  the  point  \|r*  =  1*0  in  the 
r0-plane,  then  from  figure  12 


Q 


ds  cos  P 


3u 

“  Si* 


ds  sin  p 


where 


and 


bv 

cis 


bv  dr  + 
5r  ds 


bv  d 6  bv 

3©  di-  =  S?  cos  p 


3v  .  0 

30  sin  P 


3u 


P 


(Bl) 


(Bla) 


(Bib) 


At  the  point  for  \|r*  =  1.0,  however,  u  and  v  are  equal  to  zero  and, 
because  the  streamlines  are  normal  to  the  passage  center  line, 


Su 


0 


so  that  the  continuity  equation  (l)  gives 


and  equation  (3a)  becomes 


1  bu  o  ,  bv 

F3e  =  2  +  3F 


(B2a) 


(B2b) 


(B2c) 


From  equations  (Bl)  and  (B2) 

ft  =  ^  (cos2  P  -  sin2p)  -  2  sin2p 
or,  from  equation  (2b) , 

-ft  =  A  +  2(l  -A)  sin2p  (B3) 
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Figure  2.  -  Fluid  particle 
in  re-plane. 


2617 


Direction  of  rotation 


Angle,  Q/i A0 

Figure  5.  -  Fluid-line  positions  of  central 
S2-  and  Sg-sur faces  for  example  II.  Hub- 
tip  ratio,  rh,  0.7000;  blade  spacing,  A 6, 
5°7 ’ ;  minimum  value  of  -0.00157. 
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Direction  of  rotation 

Sg-surface 


Angle,  0/A 6 

Figure  6.  -  Fluid-line  positions  of  central  Sx-  and  S2-surfaces  for  example  III.  Hub-tip 
ratio,  rh,  0.7000;  blade  spacing,  A9,  20°26,j  minimum  value  of  i,  -0,. 01318. 


Radius , 
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Direction  of  rotation 


Angle,  0/A0 


Figure  7.  -  Fluid-line  positions  of  central  and 

So-surfaces  for  example  IV.  Hub-tip  ratio,  rh, 
0.50105;  blade  spacing,  A0,  10°13 1 ;  minimum  value 
of  -0.00544. 
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Direction  of  rotation 


Z-  1  13  15  3  7 

U  "S  4^2  8  4  F 


Angle,  6/ A0 


Figure  8.  -  Fluid-line  positions  of  central  S and  S2-surfaces 
for  example  V.  Hub-tip  ratio,  r^,  0.89453;  blade  spacing,  A0, 
10o13,;  minimum  value  of  f,  -0.00231. 
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Direction  of  rotation 


Figure  9.  -  Definitions  of  terms  used  to  describe  deviation 
of  flow  surfaces  from  initial  orientation  in  re-plane  at 

a  -  1.0. 


Figure  10.  -  Variation  in  rotation  (p  -  p^  of  central  Sj-surfaces  with  angle 

Equation  (7a). 


from  their  assumed  orientation  in  two-dimensional  from  thelr  assumed  orientation  in  two-dimensional 

solutions  is  small.  solutions  is  small. 


